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The spontaneous polarization, dynamical Born charges, and piezoelectric constants of the III-V 
nitrides AIN, GaN, and InN are studied ab initio using the Berry phase approach to polarization 
in solids. The piezoelectric constants are found to be up 10 times larger than in conventional 
III-V's and II-VI's, and comparable to those of ZnO. Further properties at variance with those of 
conventional III-V compounds are the sign of the piezoelectric constants (positive as in II-VI's) and 
the very large spontaneous polarization. 
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In this paper we report an ab initio study of the sponta- 
neous polarization, piezoelectric constants, and dynami- 
cal charges of the III-V nitride semiconductors AIN, GaN, 
and InN.El This class of polarization-related properties 
is of obvious, importance for the study of nitride-based 
piezodevicesa and multilayer structures. In particular, 
the knowledge of these properties allows an insightful 
treatment of the polarization (and ensuing electric fields) 
in strained and polarized nitride junctions and superlal^j 
tices under any strain condition, as discussed elsewhere.u 
From the present study, one of the firstpapplications of the 
modern theory of polarization in solidso'cl to real and "dif- 
ficult" materials of technological interest (the first to our 
knowledge using pseudopotentials), the nitrides emerge 
as highly unusual III-V materials, resembling II- VI ox- 
ides and in some respects ferroelectric perovskites. The 
results we report here are of special interest in view of the 
scarcity of the data (both experimental and theoretical) 
available at present for the nitrides. 

In the absence of external fields, the total macroscopic 
polarization P of a solid is the sum of the spontaneous 
polarization P"'^ in the equilibrium structure, and of the 
strain- induced or piezoelectric polarization SP. In the 
linear regime, the piezoelectric polarization is related to 
the strain e by 



(1) 



which defines the components of the piezoelectric tensor 
Cij (Voigt notation is used). 

The natural structure of the III-V nitrides is wurtzite, 
an hexagonal crystal structure defined by the edge length 
a of the basal hexagon, the height c of the hexagonal 
prism, and an internal parameter u defined as the anion- 
cation bond length along the (0001) axis in units of c. 
Wurtzite is the structure with highest symmetry cot 
patible with the existence of spontaneous polarization! 
and the piezoelectric tensor of wurtzite has three non- 
vanishing independent components. Therefore, the po- 
larization in these materials system will both a sponta- 
neous and a piezoelectric component. 



In this work we shall restrict ourselves to polarizations 
along the (0001) axis; this is the direction along which 
both standard bulk materials and nitride superlattices 
are grown. The spontaneous polarization vector then 
points in the direction of the c axis, V'^'^ — P°'^z, and 
the piezoelectric polarization is simply expressed via the 
piezoelectric coefficients 633 and 613 as 



SP3 = 63363 -I- 631(61 + 62), 



(2) 



where 63 = (c — co)/co is the strain along the c axis, 
61 = 62 = (a — ao)/ao is the strain in the basal plane, ao 
and Co being the equilibrium values of the lattice parame- 
ters. The third independent component of the piezoelec- 
tric tensor, 615, is related to the polarization induced by 
a shear strain, and will not be considered in this work. 

To linear order, the change in polarization can also be 
expressed as 

0P3 ^ -^(a- ao) + ——(c~ Co) + ——{u~ uq, (3) 
oa oc ou 

where Oq, Cq, and Uq are the equilibrium lattice constants 
and internal parameter. (Note that the parameters a, c 
and u are not independent, because of the condition of 
vanishing forces along the axis direction.) Once the po- 
larization derivatives in Eq. ^ are known, the macroscopic 
piezoelectric tensor coefficients can be expressed as 



and 



where 



dPz , 4eco ^ du 

633 = Co -5 1- ^ „ Z — 

oc V3a5 dc 
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(4) 



(5) 



(6) 



is the axial component of the Born, or transverse, dy- 
namical charge tensor Z'-'^K In the last equation above. 
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it is implicit that the vector connecting the cation with 
the anion has modulus uc and points in the direction 
of the c axis. The first term in Eqs. ^ and || is called 
the clamped-ion term, and it represents the effect of the 
strain on the electronic structure. The second term quan- 
tifies the effects of internal strain on the polarization. 
The derivatives of u with respect to c and a in Eqs. ^ 
and ^ are related to the strain derivatives of u through 
Co du/dc = du/dcs and oq du/da = 2 du/dei. 

We have calculated the polarization within the Berry 
phase approach of Ref. ^ The method's output is the 
difference in polarization among two states of a system, 
provided they can be connected by an adiabatic transfor- 
mation which leaves the system insulating. The differ- 
ence of electronic polarization APe among two systems 
is then related to the geometric quantum phase by 



Pe(A2)-Pe(Ai) 



where 



(7) 



(8) 



the integration domain in momentum space is the recip- 
rocal unit cell, A characterizes the adiabatic transforma- 
tion, and t/)*-'^-' is the geometric quantum phase 



0(^)(k, k') = lm{ln[det{u(^\k)\ul'\k'))] 



(9) 



of the occupied crystal Bloch states ul^-'(k). The defini- 
tion of the geometric quantum phase introduces an arbi- 
trary constant in the value of the polarization so that the 
latter is well-defined only modulo eR/fi, with R a real- 
space lattice vector. This is not a difficulty in practice, 
as one is generally interested in polarization changes such 
that |AP| <C |eR/r2|, so that no ambiguity can arise. 

The adiabatic transformation labeled by A is trivially 
identified as far as the strains ei, £2, and £3 are concerned, 
so that the strain derivatives of the polarization within 
the wurtzite structure are easily calculated. On the other 
hand, the absolute value of the spontaneous polarization 
of the wurtzite structure (state A2 in Eq. ^ must be 
calculated as a difference relative to a structure having 
zero polarization, and functioning as reference state Ai, 
since only polarization differences can be meaningfully 
calculated within the modern theoryJ^'El: zincblende is the 
natural choice for this purpose. One needs not actually 
devise a gap-preserving adiabatic transformation turning 
zincblende into wurtzite (a quite non-trivial task): the in- 
terface theorem of Ref. g proves that Eq. ^ can be used 
directly provided that an insulating interface can be built 
connecting the_two different structures. As this is found 
to be the caseH for the systems of relevance here, the 
absolute polarization of wurtzite nitrides can be easily 
computed (see also Ref. 0). 

All the present first principles calculations are based on 
density functional theoryLl in the local density approxi- 
mation (LDA) for the exchange-correlation enjecgy func- 
tional, for which we adopt the Ceperley-Aldert3 form as 



parametrized by Perdew and Zunger.0 The wave func- 
tions are expanded in a plane-wave basis set up to-an en- 
ergy cutoff of 25 Ry. Ultrasoft pseudopotentialala have 
been employed for all the elements involved in the cal- 
culations. In particular, the Ga and In pscudopotcntials 
are constructed so as to explicitly include the semicore 
d electrons in the valence, resulting in a very accurate 
description of the electronic and structural properties. 

All calculations for the wurtzite structure used a mesh 
of 12 Chadi-Cohen special points for the Brillouin zone 
sampling.Ej The k-space integration for computing the 
polaciaation was instead done on a (4,4,10) Monkhorst- 
Packta mesh; this amounts to 16 k^ points times 10 k|| 
points. Refining the mesh to 360 points (the (6,6,10) 
mesh) gave no appreciable variation of the calculated val- 
ues. The partial derivatives of the polarization appearing 
in Eq. ^ are easily evaluated numerically through the po- 
larization change induced by typically a ±2% change in 
the crystal parameters a and c (with u kept fixed, i.e. 
clamped ion) and in the internal parameter u, around 
their equilibrium values. 

Since we are interested in the spontaneous polariza- 
tion, a careful determination of the equilibrium structure 
parameters is needed. The three independent structural 
parameters have been obtained by polynomial interpola- 
tion of the total energy values calculated over a grid of 
values of a and c. The internal parameter u was opti- 
mized at each (a,c) grid point following the Hellmann- 
Feynman forces. The results, summarized in Tablefl, are 
very close to those reported in our previous work.E3 

In Table ^ we report the spontaneous polarization, 
dynamical charges, piezoelectric coefficients (and parts 
thereof) for the III-V nitrides. Also included are the 
same quantities calculated for ZnO and BeO (with the 
Zn 3d states included in the valence). In addition, we 
list the same quantities as calculated in Ref. |l^, using 
the Berry phase approach and the FLAPW method, for 
ZnO and BeO. The agreement of our results with experi- 
ments for the latter materials is quite good. To facilitate 
further comparison with other III-V and II- VI systems, 
we have collected in Table III the piezoelectric constants 



for a number of III-V and II- VI compounds from Refs. 
ir?! , [1^ , and |l9|, pCalculated by density functional linear 
response theoryEj 

To compare zincblende and wurtzite structures in Ta- 
bles ^ and |ll, we have transformed the zincblende piezo- 
electric tensor by rotating to a coordinate system with 
z along a (111) direction. Thus, 633 = 2ei4/\/3 and 
631 — — ei4/V3, with ei4 the piezoelectric constant of 
the zincblende material referred to the cubic axes. This 
direct comparison is meaningful since (1) the possible 
deviations from the ideal wurtzite structure do not influ- 
ence significantly the polarization derivatives, and (2) it 
is known that in polytypical materials the values of the 
piezoelectric constants in the two competing structures 
agree to within a few perceat (modulo the appropriate 
coordinate transformation) Jl3 
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The III-V nitrides exhibit rather striking and unusual 
properties. Overall they resemble II- VI oxides, and ZnO 
in particular, and appear to be very different from con- 
ventional III-V semiconductors. Several features of the 
data in Tables || and III confirm this statement: 

1) The absolute value of the piezoelectric constants is up 
to about 10 times larger than in conventional III-V's and 

II- VPs. In particular, both the constants 633 and 631 
of AIN are larger than those of ZnO, and are therefore 
the largest known so far among the tetrahedrally bonded 
semiconductors . 

2) The spontaneous polarization, i.e. the polarization at 
zero strain, is very large in the nitrides. That of AIN, 
in particular, is the largest of our set of semiconductors, 
and only about 3-5|-times smaller than in typical ferro- 
electric perovskites.E]l Indeed, the piezoelectric constants 
themselves aije-|Only about 3 times smaller than in those 
ferroelectrics£3 

3) The piezoelectric constants have the same sign as in II- 
VI compounds, and opposite to III-V's. While in normal 

III- V's the clamped-ion term is larger in absolute value 
than the internal-strain ionic contribution, in the nitrides 
the latter prevails due to the larger Z* . Compared to nor- 
mal III-V's, this sign inversioi)-4s a qualitative difference 
of obvious practical relevance.E^ 

4) The nitrides follow qualitatively a well defined III-V 
trend: the piezoelectric constants increase in magnitude 
as a function of the anion chemical identity as one moves 
upwards within Period V, i.e. from Sb to N, because 
the ionic contribution tends to prevail over the electronic 
"clamped-ion" term as the anion becomes lighter (note 
that a similar trend is also obeyed by Zn II- VPs). The ni- 
trides are an extreme case of this trend, and their piezo- 
electric response is by far larger than that of all other 
III-V's (see Table [II), and of opposite sign. 

5) The Born effective charge in the nitrides is within 
10% of the nominal ionic charge, i.e. the ratio Z* /Zion 
is about unity. The same ratio is at most 0.7 in III-V 
compounds. In this context, it is appropriate to men- 
tion that the effective charges have no rigorous relation 
to nominal ionicity (although they may be a reason- 
able qualitative measure of ionicity in many binary com- 
pounds, suggesting trends analogous to those given by 
e.g. thn^ charge asymmetry g parameter of Garcia and 
Coheno) In particular, Z* needs not be smaller than 
or equal to the nominal ionic charge except in a purely 
ionic picture; conversely, whereas Z* < Z^on may be in- 
terpreted as "smaller ionicity" in some sense, Z* > Zion 
is not meaningful in an ionic picture. Effective charges 
much larger than the nominal ionic charge are common 
in perovskitesu aud (to a lesser extent) even in alka- 
line earths oxides.E£l This behavior is due to covalency 
and correlation contributions. We find (see Table ||) 
Z* > Zion for InN, but with a deviation of ~1%, i.e. the 
effective charge is again very close to nominal ionicity. 

Because of the sensitive dependence of the polarization 
on the structural parameters, there are some quantitative 
differences in spontaneous polarization for the three ni- 



trides. For instance, the increasing non-ideality of the 
structure going from GaN to InN to AIN {u gets longer 
and c/a gets shorter), corresponds to an increase of the 
spontaneous polarization. To check this, we calculated 
the spontaneous polarization of AIN using the ideal struc- 
tural parameters Uid — 0.375 and Qd = ■\/8/3a, and ob- 
tained a spontaneous polarization equal to -0.033 C/m^, 
very similar to GaN and InN, and a factor of 2.5 smaller 
than the spontaneous polarization in the actual equilib- 
rium structure. The piezoelectric constants (i.e. polar- 
ization derivatives) are instead expected to vary mildly. 

Additional Ill-V-like trends to be observed are that in 
III-V's the effective charge has similar values fo*, AlAs 
and GaAs, but a rather larger value for InAs.ll3 This 
trend is clearly obeyed by the nitrides as well, indeed in 
a more clearcut way than in most other III-V's. Another 
interesting trend is the connection between the chemical 
nature of the cation and the piezoelectric constant. The 
pieapelectric constant for Al and In III-V's was found to 
beO consistently larger than the corresponding value for 
Ga compounds. This is also true of the nitrides, where 
the deviations are again larger. Since our calculations 
include the semicore d states in the valence for both In 
and Ga, this trend is not an artifact of the neglect of 
semicore states as hypothesized in Ref. |l^. 

While it is known that treating the semicore d states as 
valence is necessaoz: ±0 obtain good structural properties 
for GaN and InN,tj'L3 neither the spontaneous polariza- 
tion nor the effective charges are influenced strongly. Us- 
ing a Ga pseudopotential in which the 3d states had been 
frozen in the core, we recalculated the spontaneous polar- 
ization and effective charges of GaN using the structural 
parameters given in Table |: the results are P°'^ ~ —0.031 
C/m^ and Z* — 2.73, showing that the direct influence 
of the 3d states on the final polarization value is minor. 
This was to be expected since the d bands in GaN are 
quite flat over most of the Brillouin zone, so that the 
overlap matrix (Eq. ^ is nearly the identity for those 
states, and its contribution to the Berry phase is small. 

Experimental data to compare our predictions with 
are very scarce. No data is available for the sponta- 
neous polarization, a quite elusive quantity in this re- 
spect; our calculated values will be a useful input for 
device modeling.cHa Also in the case of the piezoelec- 
tric constants, our predicted data should be of use to 
future ejfp(erimental work. The only experimental data 
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availablcEZl to our knowledge are for AIN, 



-^33 



C/m^, and Cg^^ = —0.58 C/m^. Our predicted values 
of 1.46 C/m^ and -0.60 C/m^ agree with experiment to 
within ~ 5%. 

In summary, we have calculated the spontaneous po- 
larization, dynamical charges, and piezoelectric constants 
of III-V nitrides by means of the Berry phase approach. 
The results show that III-V nitrides resemble II- VI com- 
pounds in terms of the sign of the spontaneous polariza- 
tion and of the piezoelectric constants; the latter con- 
stants are much larger in absolute value than those of all 
III-V's and II- VPs, except for ZnO. The values of sponta- 
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neous polarization are comparable to or larger than those 
of II-VI's; AIN has the largest value reported so far for 
any binary compound, and only a factor of 3 to 4 smaller 
than in typical perovskites. Some remnant of normal 
Ill-V-like behavior survives, embodied in the trends of 
the piezoelectric constants and effective charges when the 
chemical identity of the anion or cation changes. 
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TABLE I. Structural parameters for AIN, GaN and InN. 





ao (bohr) 


Co /ao 


uo 


AIN 


5.814 


1.6190 


0.380 


GaN 


6.040 


1.6336 


0.376 


InN 


6.660 


1.6270 


0.377 



TABLE II. Calculated spontaneous polarization (in units 
of C/m'^), Born effective charges (in units of e), and piezoelec- 
tric constants (in units of C/m^) for III-V wurtzite nitrides 
and II- VI wurtzite oxides. 





pcq 


Z* 


du/des 


633 


631 


633 


631 


AIN 


-0.081 


-2.70 


-0.18 


1.46 


-0.60 


-0.47 


0.36 


GaN 


-0.029 


-2.72 


-0.16 


0.73 


-0.49 


-0.84 


0.45 


InN 


-0.032 


-3.02 


-0.20 


0.97 


-0.57 


-0.88 


0.45 


ZnO 


-0.057 


-2.11 


-0.21 


0.89 


-0.51 


-0.66 


0.38 


BeO 


-0.045 


-1.85 


-0.06 


-0.02 


-0.02 


-0.60 


0.35 


ZnO" 
BeO" 


-0.05 
-0.05 


-2.05 
-1.72 


-0.25 
-0.09 


1.21 
0.50 


-0.51 


-0.58 
-0.29 


0.37 
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TABLE III. Piezoelectric constants (in units of C/rr?) of 
several zincblende compounds, calculated via linear response 
theory in Refs. (IH-V), |l| (CdTe) and |l9| (other ILVI). 





633 


631 




633 


631 


CdTe 


0.03 


-0.01 


AlAs 


-0.01 


0.01 


ZnS 


0.10 


-0.05 


GaAs 


-0.12 


0.06 


ZnSe 


0.04 


-0.02 


InAs 


-0.03 


0.01 


AlP 


0.04 


-0.02 


AlSb 


-0.04 


0.02 


GaP 


-0.07 


0.03 


GaSb 


-0.12 


0.06 


InP 


0.04 


-0.02 


InSb 


-0.06 


0.03 
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